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ABSTRACT 

We analyse the kinematics and chemistry of the bulge stars of two simulated disc 
galaxies using our chemodynamical galaxy evolution code GCD+. First we compare 
stars that are born inside the galaxy with those that are born outside the galaxy and 
are accreted into the centre of the galaxy. Stars that originate outside of the bulge 
are accreted into it early in its formation within 3 Gyrs so that these stars have high 
[a/Fe] as well as having a high total energy reflecting their accretion to the centre of the 
galaxy. Therefore, higher total energy is a good indicator for finding accreted stars. The 
bulges of the simulated galaxies formed through multiple mergers separated by about 
a Gyr. Since [a/Fe] is sensitive to the first few Gyrs of star formation history, stars 
that formed during mergers at different epochs show different [a/Fe]. We show that the 
[Mg/Fe] against star formation time relation can be very useful to identify a multiple 
merger bulge formation scenario, provided there is sufficiently good age information 
available. Our simulations also show that stars formed during one of the merger events 
retain a systematically prograde rotation at the final time. This demonstrates that the 
orbit of the ancient merger that helped to form the bulge could still remain in the 
kinematics of bulge stars. 

Key vifords: Galaxy: bulge — Galaxy: kinematics and dynamics — galaxies: Inter- 
actions — galaxies: Formation — galaxies: evolution 



1 INTRODUCTION 

In recent years, there has 
disc galaxy formation in a 



been progress in modelling 
Cold Dark Matter (CDM) 
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; Brook et al.. ,2004: Governato et al.l 


20071: lOkamoto et al. 


2008|). Some studies include both Type 



II (SNe II) and Type la (SNe la) supernovae and discuss the 
details of the chemical pro perties in the simu lated gala xies in 
an isolated halo collapse (iRaiteri et al.l 1996: Bcrczi^| l999l) 



or hi erarchical clustering l|Brook et al.l l2005i : iRenda et al.l 
I2OO5I ). These studies mainly focus on t he disc component of 
the si mulated galaxies. None, except iNakasato fc Nomotd 
l|2003l ). study the bulge for mation in three dime n siona l 
chemodynamical simulations. INakasato fc Nomotd (120031 ) 
suggest that the bulge may be composed of two chemically 
different components, one being formed through an early 
epoch subgalactic clump merger, and the other in the in- 
ner region of the disc after this merger. This paper focuses 
on the population of bulge stars built up through mergers. 
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using higher resolution and more sophisticated chemody- 
namics simulations. Especially, we study the properties of 
accreted and locally forme d stars within the bulge, which 
INakasato fc Nomotd l|2003l ) do not look at. 

There are two major formation scenarios for bulge stars. 
Th e first scenario is secular evolution of the disc (see , 
e.g. ICombes et al]|l990l : iNorman et al]|l996l : lNoguchi|[l998l ) . 
In most secular evolution scenarios, gas and stars from 
the disc are tunneled to the central regions of the galaxy. 
ICombes et all (|l990l ) showed that it is possible to form a 
bulge via the secular formation of a bar, and iNoguchi (, 19981 ) 
showed that the disc of a young galaxy may become grav- 
itationally unstable and fragment into large pieces. These 
large clumps fall into the central region due to dynamical 
friction a nd proceed to form a m a ssive bulge via mergers. 
iNoguchil (|T99S ) and llmmeli et al] l|2004l ) also suggest that 
several peculiar morphological structures seen in the Hub- 
ble Deep Field can be explained by a fragmented galactic 
disc model. 

The other major formation scenari o for bulge stars is 
via m ergers of sub clumps, first shown bv lToomre fc Toomrd 
1 19721 ). Using numerical simulations and i ncluding star 
formation, gas dynamics and SNe feedback, iKat j (|l992l ) 
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Table 1. Simulation parameters 



Name 








niDM 


^gas 


EDM 




ho 






(Mq) 


ikpc) 


(Mq) 


{Mq) 


ikpc) 


(kpc) 








Gall 


8.8 X 10" 


240 


9.2 X IQS 


6.2 X 10^ 


0.57 


1.1 


0.3 


0.7 


0.039 


Gal2 


9.3 X 10" 


270 


3.3 X 10*^ 


1.9 X 10^ 


0.87 


1.5 


0.3 


0.65 


0.045 



showed that hierarchical clustering predicts a number of 
mergers at high redshift. In these merger processes, the num- 
ber of bulge stars builds up by accreted stars which were 
already formed within the infalling galaxies, and by stars 
which formed locally within the main galaxy during and 
after the merger via starbur sts. Although a number of stud- 
ies focus on accreted s tars (|Brook et al.l l2003l : iHelmi et al.l 
l2006l : iFont et al.|[200^ '). there have been few studies com- 
paring the properties of these accreted satellite stars with 
those of the s tars f ormed locally in the galaxy. Recently 
IZolotov et al.l (|2009l ) compared the properties of accreted 
and locally formed stars in the halo of the galaxy. There are 
very few studies however, dedicated to the bulge chemody- 
namical properties. This is the first study to compare ac- 
creted and locally formed stars within the bulge. 

The chemical properties of stars are important to under- 
stand the formation history of the bulge. Chemical elements 
heavier than He are the end products of stellar evolution. 
The so-called a-elements and iron (Fe) are of particular 
interest since it is known that they are produced primar- 
ily in SNe II and SNe la respectively. SNe la and SNe II 
have different timescales and thus studying the abundance 
ratios of the a elements with respect to Fe gives unique 
fossil information on the past conditions and evolution of 
the gal axy. Pure chemical evolutio n studies of the bulge 
such as iMatteucci fc Brocatd l|l99(t) : iFerreras etahl (|2003l ') 
are useful for reconstructing the star formation history of 
the bulge. IMatteucci fc Brocatd (|l990h studied the chemi- 
cal abundances in the bulge, and looked at the [a/Fe] ratios 
for different elements. They predicted a relative enhance- 
ment of [a/Fe] in the bulge, due to the high Star Formation 
Rate (SFR). Based on the observed distributi on of metallici- 
ties an d using a simple star formation model, IFerreras et al.l 
l|2003h found that very short infall timescales are required 
in the bulge. These studies are useful to reconstruct the 
star formation history of the bulge. However, pure chemical 
evolution models do not answer how such a star formation 
history occurred in the galaxy, or predict any information 
about kinematics. 

The findings of IMatteucci fc Brocato] lll99(t) were con- 
firmed by the observat io nal s tudies of Mc William fc RichI 
||1994 ). iFulbright et~al] (120061 ) . and IZoccah et al.l ll2006f). 
By lo oking through Baade's Window, Mc William fc Rid] 
l| 19941 ) found that the primary enrichment process of the 
bulge was SNe II. Recently, 8m class large telescopes have 
enabled observers to observe stars in the Milky Way bulge 
in mor e detail. Using a larg er sample of stars at higher res- 
olution IZoccah et aLl (|2006') concluded that the Milky Way 
bulge likely formed through a short series of starbursts trig- 
gered by gas rich mergers in the early universe. In the future, 
both kinematical and chemical properties of bulge stars will 
become available via projects such as APOGEE which can 
look through the obscuring dust to the bulge region of the 



Milk y Way using n ear infrared multi-object spectrographs 
(AUende Prieto et al..2008 ). ESA's Gaia mission should pro- 
vide six dimensional phase space and chemical information 
on up to one billion stars in the Milky Way. These stars 
lie predominantly in the disc however there is some hope 
that G aia may also ob serve significant numbers of bulge 
stars (jRobin et aLlbOOSl ). Uhimately, JASMINE, a planned 
Japanese space astrometry mission, will observe around ten 
million stars in the bulge b y taking detailed astrometric 
readings in the infrared band (|Gouda et al.|[2006l ) . This huge 
influx of new data is expected to spur new scientific results in 
different fields of astronomy. However, to extract the forma- 
tion history of galaxies from such a large amount of observa- 
tional data, theoretical models and predictions are crucial. 
Therefore, we initiate a project to study how we can tell 
the past formation history of the bulge through the current 
properties of bulge stars. 

We study both chemical and kinematical properties 
of bulge stars using ACDM cosmological simulations. We 
take two high res olution simulations from the sample of 
iBailin et al.l |2005h . In our ACDM simulation, bulges formed 
during hierarchical clustering of building blocks at high red- 
shift. We therefore study the kinds of properties that are 
expected in bulge stars, if the bulge was formed through 
such a hierarchical clustering scenario. Therefore, our study 
does not mean to reproduce the Milky Way bulge which may 
or may not form through hierarchical clustering. 

The outline for this paper is as follows. In Section 2 
we describe our numerical simulations and define our bulge 
stars. In Section 3 we present the results of our chemody- 
namical analyses on accreted and locally formed stars within 
the bulge. Finally we present our conclusions in Section 4. 



2 THE CODE AND MODEL 

To simulate our galaxies, we use our original galac- 
tic chemodynami cs evolution code GCD+ developed by 
iKawata fc GibsonI (|2003. ). GCD+ is a three-dimensional tree 
A^-bo d y /smoothed particle hydro dynamics (SPH 



19771; [Cingold fc MonaghanI 



19861 : iHernguist fc Katzlll989l : 



Lucy 



19771) code (iBarnes fc Hut 



Katz et al]ll996l ) that incor- 



porates self-gravity, hydrodynamics, radiative cooling, star 
formation, supernova feedback, and metal enrichment. GCD+ 
takes into account chem ical enrichm ent by both SNe II 
JWooslev fc WeayeJll995l ) and SNe la (|lwamoto et al.lll999l : 
iKobavashi et all I2OOOI) and mass loss from intermediate- 
mass stars l|van den Hoek fc GroenewegenI ^1997^ . and fol- 
lows the chemical enrichment history of both the stellar and 
gas components of the system. 

Radiative cooling, which depend s on the metallicity of 
the g as (derived with MAPPINGSIII: ISutherland fc Dopital 
Il993h is taken into account. The cooling rate for a gas with 
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Figure 1. [Fe/H] (upper) and [Mg/Fe] (lower) against formation 
time (tf) for Gall (left) and Gal2 (right). Red filled and blue 
open circles are in situ and accreted stars respectively. There are 
3 major periods of star formation clearly distinguishable in Gall. 
Star formation is smoother in Gal2. 



solar metallicity is larger than that for gas of primordial 
composition by more than an order of magnitude. Thus, 
cooling by metals should not be ignored in numerical sim- 
ulations of gala xy formation (iKallander HultmanI flQQSl : 
iKav et al]l2000l ). However, we ignore the effect of the UV 
background radiation for simplicity. 

The two galaxie s simulated here are from the sample 
of lBailin et alTl|2005l ), KGCD and AGCD, hereafter "Gall" 
and "Gal2" respe ctively. Gall i s a hi gh resolution version 
of galaxy "Dl" in lKawata et al.l (|2004l '). We used the multi- 
resolution technique in order to maximise the mass resolu- 
tion wi thin the regions whe re the disc progenitors form and 
evolve (jKawata et al.|[2004l ). The initial conditions for Gal2 
were kindly provide d by M. Steinmetz, and are described in 
lAbadi et al.l l|2003al lbh. 

We summari se the pro perties of the galaxies in Table 
1 taken from Connors et al.l (200^). The first column repre- 
sents the galaxy name; the second column, the virial mass; 
Column 3, the virial radius; Columns 4 and 5 represent the 
mass of each gas and DM particle in the highest resolution 
region, and Columns 6 and 7 are the softening lengths in that 
region. The cosmological parameters for the simulation are 
presented in Columns 8-10. flo is the total matter density 
fraction, ho is the Hubble constant (100 kms~^Mpc~^) and 
fib is the baryon density fraction in the universe. Note that 
the cosmology is slightly different between the two models, 
and the age of the Universe is 13.5 and 14.5 Gyrs for Gal 1 
and Gal2 models respectively. 

To identify the main progenitor galaxy, at each output 
of the simulation, a friend-of-friend (FOF) group finder is 
used to identify stellar groups. Specifying a linking length 
b and identifying all pairs of particles with a separation 
equal to or less than b times the mean particle separation 
as friends, stellar groups are defined as sets of particles con- 
nected by one or more friendship relations. The other main 
parameter in FOF is the minimum number of particles. By 
setting this parameter sufficiently high, one avoids including 
spurious objects that may arise due to chance. In our simu- 
lations we use a linking parameter b = 0.01 and a threshold 
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Figure 2. Global star formation rate (SFR) for Gall (left) and 
Gal2 (right). Red solid, blue dashed and black dotted are in situ, 
accreted and total stars respectively. 



number of particles of 100. We define the largest group which 
has the highest number of the FOF identified particles as the 
main progenitor galaxy. 

Note that the vers ion of GCD+ used in this paper 
llKawata fc Gibson 20031) applies the SNe H yields calcu- 
lated bv lWooslev fc Weaver! (|l995l ). It is well known that this 
yield seems to be overestimated and leads to lower [a/Fe], 
compared with that obs erved in low-metal l icity stars in 
the solar neighbour hood (|Timmes et al. I I 19951 : lGibsonlll99'?1 : 
I Gibson et al.l Il997l ) . We adopted pure Woosley & Weaver 
yields for Gall. On the other hand, we apply a half Fe yield 
for Gal2. As a result, [Mg/Fe] values in Gall are underesti- 
mated compared to Gal2. Therefore, we cannot compare the 
absolute value of [Mg/Fe] between Gall and Gal2. However, 
we compare the relative difference of [Mg/Fe] among dif- 
ferent samples of stars within a simulated galaxy. Our code 
does not "disperse" metals between gas particles. Therefore, 
the spread in the metallicity distributions will be artificially 
high, while the peak of the distributions (the average) should 
be robust (or at least the distances between the peaks) . 

We identify bulge particles using the simulation output 
at the latest time. For Gall, we use the output at 2 = 0.1, 
as going to any lower redshift results in an unacceptable 
amount of contamination from low-resolution particles in 
our simulated galaxy. For Gal2, we use the output at 2 = 0. 
To define the bulge, we first set the disc plane of our galax- 
ies to be the x-y plane and the rotation axis to be the 
z-axis. The bulge was defined as the central 2 kpc of the 
galaxy, excluding ±1 kpc in the z-direction due to contami- 
nation with disc stars. Avoiding the disc plane will provide 
us with a "cleaner" sample since the central region of the 
disc suffers from "overcooling" problems and therefore pro- 
duces too many stars throughout the simulations. All stars 
found within the remaiinng volume described above we refer 
to as "bulge stars", and we analyse their properties in the 
next section. 



3 RESULTS 

Fig.[T]shows [Fe/H] (upper panel) and [Mg/Fe] (lower panel) 
as a function of the time (tf) when bulge stars in our two 
galaxies formed. There is clearly a different formation his- 
tory for the bulges of the two galaxies. The bulge of Gall 
forms via a series of three mergers within the first 5 Gyrs. 
The bulge of Gal2 is formed in a similar timescale, however 
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Figure 3. [Mg/Fe] vs [Fe/H] for Gall{left) and Gal2 (right). 
Red filled and blue open circles are in situ and accreted stars 
respectively. As expected from Fig.[T]accreted stars have generally 
higher [Mg/Fe] ratios. A population of accreted stars in Gal2 have 
high [Mg/Fe] because these stars formed in a large satellite which 
merged into the main galaxy at an early epoch. 
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Figure 4. Total energy (i?tot) vs angular momentum (iz) for 
Gall{left) and Gal2 (right) for accreted (open blue circles) and 
in situ (filled red circles) stars. Accreted stars have higher Etot- 
For Gal2 the distinction is unclear. This is due to a major merger 
early in its history bringing stars of a satellite galaxy into the 
main galaxy as accreted stars. 



without the distinctly separated star formation epochs. Star 
formation here is more smooth and continuous. Both galax- 
ies formed through hierarchical clustering, and bulge forma- 
tion is associated with a series of mergers at early epochs. 
However, Gall and Gal2 have different merger histories and 
therefore star formation histories. The aim of this study is to 
find out how we can tell such a difference in merger histories 
at a high redshift from the current properties of bulge stars. 
To this end, first we separate the bulge stars into accreted 
stars and locally formed stars, and compare their properties 
in Section 13.11 From the current chemical and kinematical 
properties of bulge stars, we study how we can distinguish 
accreted stars from locally formed stars. In Section 13.21 we 
compare the properties of the bulge stars with different for- 
mation epochs to see whether the memory of mergers at 
different epochs persists or not. 

3.1 Accreted and in situ stars 

In this section, we compare the chemical and kinematical 
properties of accreted and locally formed stars. We trace 
back the formation time and location for all the bulge stars, 
and any stars that are formed locally within a radius of 5 
kpc from the largest progenitor are given the title "in situ 
stars". Stars that form at a radius greater than 5 kpc from 
the centre of the galaxy and end up within the area defined 
above as the bulge, we term "accreted stars". We chose 5 
kpc arbitrarily. We experimented using a larger radius up to 
20 kpc, and generally found the same conclusion. 

Fig- U shows [Fe/H] and [Mg/Fe] vs formation time for 
accreted and in situ stars with different symbols. Mg is one 
of the so-called Q-elements. These elements are primarily 
produced when massive stars with short lifetimes explode 
as SNe II. Fe is produced predominantly in SNe la, from 
lower mass binary stars with longer lifetimes. Accreted stars 
have lower [Fe/H] and higher [a/Fe] since they form in early 
epochs before the enrichment from SNe la becomes impor- 
tant. In situ stars continue to undergo strong star formation 
for up to 5 Gyrs (as seen in Fig. [2j which shows the SFR 
history for in situ and accreted stars). Fig. [3] shows [Mg/Fe] 
vs [Fe/H]. Note that there exists a distinct population of ac- 
creted stars in Gal2 with high [a/Fe]. This is because these 



accreted stars were formed in a large companion galaxy at 
early epochs {t < 1.5 Gyrs) that underwent a major merger 
with the main galaxy (Fig. [2} . This demonstrates that the 
probability distribution of high [Mg/Fe] stars may be sensi- 
tive to merger history at high redshift. 

We find that Fig. |4] is the most useful for identifying 
the accreted stars. It shows the stars total energy (iStot) 
plotted against their z-component of angular momentum 
(Lz) which is also used to identify acc r eted stellar groups 
in halo stars (iHelmi fc de Zeeuwl I2OO0I : iKnebe et alTl2005l : 
iFont et "al]|2006l ). The accreted stars have a relatively higher 
Etot and are spread to a large range of Lz . This is especially 
the case for Gall. For Gal2 however, the distinction is less 
clear. There exists a significant population of accreted stars 
with low Etot in Gal2. This is due to the major merger this 
galaxy experienced at an early epoch in its history. These 
accreted stars with low Etot were actually formed in the 
centre of a large satellite galaxy that merged with the main 
galaxy early on which explains why their Etot is comparable 
to the in situ stars. 

We conclude that the i?tot - Lz diagram is a powerful 
tool to identify accreted stars in the bulge. Its use is known 
for halo stars (Brook et al 2003, Morrison et al 2009). Our 
study shows that it can be applied for bulge stars which 
experienced more violent merger histories. However there 
are significant overlaps with in situ stars. Especially, it is 
difficult to identify stars which are born in a large satellite 
galaxy that merged at an early epoch. 



3.2 Multiple merger bulge formation 

The left panel of Fig. [2] clearly shows that Gall has 3 bursts 
of star formation. From snapshots, we confirmed that each 
of the bursts corresponds to a merger. We therefore divided 
the bulge stars into 3 samples corresponding to each of the 
bursts. The stars formed during the period of the first burst 
{tf < 1.5 Gyrs) are referred to as "old" bulge stars; those 
formed during the second burst (1.5 < i/ < 3 Gyrs) are 
called "intermediate" bulge stars; and those formed during 
and after the third and largest burst {tf > 3 Gyrs) are called 
"young" bulge stars. From the right panel of Fig. [2] we see 
that Gal2 also has three bursts of star formation, although 
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Figure 5. As in Fig.|4l but now using different aged bulge stars 
for Gall (left) and Gal2 (right). Red open circles, blue filled cir- 
cles and black crosses are old, intermediate and young bulge stars 
respectively. Trends are visible, especially in Gall which contains 
an intermediate aged population with a distinctly high L^, pos- 
sibly indicating the prograde orbit of a merger at that epoch. 
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Figure 6. As in Fig. [3] but now using different aged bulge stars. 
Red open circles, blue filled circles and black crosses are old, in- 
termediate and young bulge stars respectively. For Gall again, 
the distinction between the different groups is clearer. 



there is significant star formation between tlie bursts. Tliere- 
fore, we also assign 3 periods of star formation for Gal2. Our 
definition for the 3 time periods here was slightly different to 
Gall, reflecting the different star formation history in Gal2. 
The formation periods of old, intermediate and young stars 
for Gal2 are defined as tf < 2 Gyrs, 2 < tf < 4 Gyrs, and 
t/ > 4 Gyrs respectively. 

In Fig. [5] we replot Etot vs for the three different 
groups of bulge stars. In Gall (left panel), the old stars are 
generally clustered about the centre with low |Lz| and iJtot- 
Young stars have higher Etot and are less clustered about 
low values of \Lz\. In Gall we clearly see that the major- 
ity of intermediate stars have high L^. This indicated that 
the intermediate population may have resulted from a spe- 
cific merger and are orbiting in a prograde direction. This 
demonstrates that some information from the ancient his- 
tory of mergers can persist in Etot and properties. If 
then, in addition, the ages of the stars are measured very 
accurately, the diagram of Etot and may become a pow- 
erful tool to reconstruct the merger history. The picture for 
Gal2, again is less clear. The three populations in Gal2 are 
well spread on the diagram. No group occupies a distinct 
region. 

Fig.|6]shows that the old stars have the highest [Mg/Fe]. 
The younger stars have lower [Mg/Fe] and a higher metal- 
licity as expected. This is because at later times the con- 
tribution from SNe la has become significant. The pattern 
is less clear for Gal2, reflecting its more continuous forma- 
tion history, which diffuses out any trends. We plot the 3 
samples of bulge stars as a probability distribution function 
for [Mg/Fe] as in Fig[71 We see that they each occupy a 
separate distinct region, although there is some overlap. As 
mentioned in Section 2, our simulations do not take into 
account metal mixing between the gas particles. We overes- 
timate the scatter of metal abundance and abundance ratios 
for star particles. Therefore, if the correct metal mixing is 
taken into account, [Mg/Fe] distributions of these 3 groups 
of stars will be more clearly different. Note however that 
here we have assumed that we have accurate age informa- 
tion available for bulge stars. However, in reality, with the 
current facilities, it is difficult to obtain the ages of the bulge 
stars. Our study demonstrates that if the accurate age infor- 
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Figure 7. Histogram of [Mg/Fe] for different age samples of 
bulge stars for Gall (left) and Gal2 (right). Red solid, blue dashed 
and black dotted are old, intermediate and young bulge stars 
respectively. 

mation becomes available (as it will in the future), we can 
identify bulges formed through distinct multiple mergers us- 
ing the age-[Mg/Fe] relation (Fig. [l|. Then, for bulge stars 
that are formed as a result of mergers at different epochs, we 
would expect differing values for [Mg/Fe], and for younger 
stars to have lower [Mg/Fe]. In addition, the Etot - dia- 
gram can give us some information on the orbits of merging 
building blocks. 



4 SUMMARY 

In this study we have analysed the chemistry and the dy- 
namics of the bulge stars of two simulated disc galaxies. The 
two galaxies had different formation histories and were sim- 
ulated under different cosmologies. However both galaxies 
were similar in mass and size to the Milky W ay, and con- 
tained distinct gas and ste llar disc components (jBailin et al.l 
l2005l : IConnors et aLllioO^ . 

We started by defining the bulge region to be the cen- 
tral 2 kpc of the galaxy excluding ±1 kpc in the z-direction 
to avoid contamination with disc stars. In this study we 
compared the properties of accreted and in situ stars within 
the bulge. Since accreted stars tend to form in early epochs, 
they have lower [Fe/H] and higher [Mg/Fe] ratios. We also 
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found that even within the bulge, the accreted stars tend to 
have higher iJtot although there is significant overlap with 
in situ stars. We found that accreted stars formed in a large 
satellite galaxy that merged into the main galaxy can have 
low Etot and are difficult to distinguish from in situ stars. 

The bulge of one of the simulated galaxies was formed 
via a series of three distinct mergers in the first 5 Gyrs as 
confirmed by snapshots of the simulations. We found that 
such multiple star formation bursts will show up in the age- 
metallicity relation or more evidently in the age-[Mg/Fe] 
relation, provided there is accurate age information avail- 
able. Then, stars that formed in earlier mergers tend to 
have higher [Mg/Fe] , and different mergers have distinctively 
different [Mg/Fe]. Therefore, the probability distribution of 
[Mg/Fe] is sensitive to the number of mergers and the epochs 
of the mergers. We found that even in the bulge Etot and 
information persists reasonably well till 2 = 0. Therefore, if 
we could pick up stars formed in a merger at a specific epoch 
using the age information, Etot and would be able to tell 
us about the orbital information of the merger, like whether 
it was in a prograde or retrograde direction. 

The other simulated galaxy had undergone a major 
merger at a very early epoch, followed by a series of mi- 
nor mergers between two later major mergers, which led to 
a smoother star formation history. This meant that much 
of the chemodynamical information had been smeared out. 
We found that a significant number of accreted bulge stars 
had high values of [Mg/Fe] , due to significant star formation 
in two large progenitors at early epochs. This sensitivity to 
the star formation rate at early epochs means that finding 
such bulge stars with high values for [Mg/Fe] could depend 
on the size of the progenitors at the very early epochs. 

Despite differences, our two simulated galaxies share 
some common formation features. There is evidence for the 
bulge of both galaxies forming via three major bursts of star 
formation. This is because in hierarchical clustering scenar- 
ios, bulges are likely to form through multiple rather than 
single mergers. If then, such merger epochs are separated by 
enough time, we may obtain information on these mergers 
from the properties observed in our simulated galaxies at 

2 = 0. 

It is not yet possible from the current observations to 
obtain accurate age information for bulge stars. However 
in the near future, when increasingly accurate age informa- 
tion becomes available from combinations of the near in- 
frared astrometry satellite JASMINE, and high resolution 



spectroscopy in the optical |Fulbright et al.ll2006D and near 
infrared with APOGEE (see AUcndc Prieto et al]|2008t ). our 
findings will prove useful for identifying a possible scenario 
for the past formation history of the Galactic bulge. 
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